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We study the proton magic O, Ca, Ni, Zr, Sn, and Pb isotope chains 


from the proton drip line to the neutron drip line with the relativistic con- 
tinuum Hartree-Bogoliubov (RCHB) theory. Particulary, we study in detail 
the properties of even-even Ca isotopes due to the appearance of giant ha- 
los in neutron rich Ca nuclei near the neutron drip line. The RCHB theory 
is able to reproduce the experimental binding energies E, and two neutron 
separation energies Son very well. The predicted neutron drip line nuclei are 
280, Ca, SNI, 6Zr, 16Sn, and 2°°Pb, respectively. Halo and giant halo 
properties predicted in Ca isotopes with A > 60 are investigated in detail 
from the analysis of two neutron separation energies, nucleon density distri- 
butions, single particle energy levels, the occupation probabilities of energy 
levels including continuum states. The spin-orbit splitting and the diffuseness 
of nuclear potential in these Ca isotopes are studied also. Furthermore, we 


study the neighboring lighter isotopes in the drip line Ca region and find some 
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possibility of giant halo nuclei in the Ne-Na-Mg drip line nuclei. 
PACS: 21.10.Gv, 21.60.-n, 21.60.Jz 
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I. INTRODUCTION 


The study of exotic nuclei, which are so-called due to their large N/Z ratios (isospin) and 
their interesting properties, e.g. halo and skin, has attracted world wide attention [1,2]. With 
the recent developments in accelerator technology and detection techniques, it has come into 

_ reality to produce these exotic nuclei and study their detailed properties with the radioactive 
N a ion beam (RIB) facilities. Since Tanihata et al. discovered the first case of halo in an exotic 


~ nucleus !Li with RIB in 1985 [3], more and more exotic nuclei have been investigated with 


' various modern experimental methods to understand this attractive phenomenon better. 
) For nuclei far from the 6-stability valley and with small nucleon separation energy, the 
— valence nucleons in exotic nuclei extend over quite a wide space to form low density nuclear 
Y , matter. It is expected that the ”halo” in exotic nuclei provides an interesting case to study 
re the nuclear environment in astrophysics in the laboratory. 

Very neutron-rich nuclei and, in particular, those near the neutron-drip lines and near 


2 closed shells, play an important role in nuclear astrophysics. Their properties such as binding 


il 


© energies, neutron separation energies, deformation parameters, etc., strongly affect the way 
~ neutron-rich stable isotopes are formed in nature by the so-called r process. 

Furthermore, the exotic nuclei are expected to exhibit some other interesting phenomena 
such as the disappearance of traditional shell gaps and the occurrence of new shell gaps, 
which result in new magic numbers. Ref. [4] has reported that the new magic number N = 16 
appears in the light neutron drip line region. Ozawa et al. have suggested the mechanism 
for the formation of the new magic number, which is intimately related to the neutron halo 
formation and is common in nuclei near neutron drip line. 

It is very helpful to use self-consistent microscopic model to study the properties of exotic 


nuclei. During the past two decades, the relativistic mean field (RMF) theory has received 


wide attention, because the RMF theory is very successful in describing many nuclear phe- 
nomena for nuclei even far from stability as well as the stable nuclei. Compared with the 
non-relativistic mean field theory, RMF can reproduce the right nuclear incompressibility 
coefficient and saturation properties (Coester line) in nuclear matter and gives naturally the 
spin-orbit coupling potential. The reviews on RMF theory are given in Ref. [5-7] . The 
starting point of the RMF theory is the Lagrangian which describes the nucleons as Dirac 
spinors moving in a mean field, composed of the interaction between nucleons (protons and 
neutrons) and mesons (c,w, p), and also Coulomb field. From this viewpoint, the RMF 
theory is more microscopic in a sense of describing the nuclear system at the meson level. 
Usually for exotic nuclei, their Fermi surfaces are very close to the continuum thresh- 
_ old. In these cases, the valence nucleons could be easily scattered to the continuum states 
< due to the pairing correlation. Thus, theories which can properly handle the pairing and 


N continuum states are needed to describe the properties of exotic nuclei. For the pairing 


) interaction, the simple BCS method and general Bogoliubov quasi-particle transformation 
» are two candidates. The former is very useful and successful for stable nuclei. However, 
—< when it is extended to exotic nuclei, the occupation probability would become finite for 


ZL those continuum states and involve unphysical fermion gas In this case, we should use the 


i Bogoliubov transformation to handle the pairing correlation in exotic nuclei instead of the 

= imple BCS method. Based on the RMF theory, the relativistic Hatree-Bogoliubov (RHB) 

5 equation can be derived by quantization of the meson fields as well as the nucleon fields 

: in Lagrangian density [8]. Furthermore, in order to describe self-consistently both the con- 

© tinuum and bound states and the coupling between them, the RHB theory must be solved 
in coordinate space, i.e., the newly developed Relativistic Continuum Hartree-Bogoliubov 
(RCHB) theory [9,10]. 

The RCHB theory was extremely successful in describing the ground state properties of 
nuclei both near and far from the {-stability line. A remarkable success of the RCHB theory 
is the new interpretation of the halo in "Li [9] and the prediction of the exotic phenomenon 
as giant halos in Zr (A > 122) isotopes [11]. The giant halos are very interesting phenomena 
in exotic nuclei because the halos are formed by more than two neutrons scattered as Cooper 


pairs to the continuum region. However, the exotic Zr (Z = 40) isotopes with N > 82 are 


rather heavy for observation by the RIB facilities at the present. With the present RIB 
techniques, light drip line nuclei are more probable to be accessed with the transfer reaction 
or other reaction mechanism. It is thus extremely valuable for us to investigate the giant 
halo phenomena in lighter nuclei like Ca (Z = 20) isotopes [12]. 

The Calcium isotope chains have received much attention due to its rich experimental 
results on binding energy, density distribution, single particle energy, radius, etc. Though 
those data are now limited near the stability line, it is useful to calculate these quantities 
in order to test the microscopic theory by future experiments. Also, there lie shell effects 
in this chain due to the short shell period as the traditional magic numbers are N = 20, 
28, 50 and sub-magic number N = 40. The properties of Ca isotope have been investigated 
with various methods based on the mean field theory, e.g., the non-relativistic Hartree- 

> Fock (HF) and Hartree- Fock Bogoliubov (HFB) method [13], the Skyrme Hartree Fock 
N (SHF) method [14], relativistic density dependent Hartree-Fock method [15], and the self- 


consistent Hartree-Fock calculation plus the random-phase approximation [16]. We apply 
oð here the RCHB theory to investigate the ground state properties of the Ca isotope chain, 
<= especially probing the halo properties in the exotic nuclei near the neutron drip line. 


In the former letter [12], we have reported briefly the halos discovered in Ca isotopes 


` = near the drip line region with the RCHB theory. Here we give details of the investigation 
= of the mechanism of the appearance of the giant halos. Besides the detailed discussion on 
3 Ca isotopes, the ground state properties of other proton magic isotope chains, O, Ni, Zr, 
"= Sn, and Pb are discussed in this article. The paper is organized as follows: In Sec. II, a 
O brief outline of the RCHB theory is given. In Sec. III, we provide the numerical results on 
these proton magic even-even nuclei and the miscellaneous properties for Ca isotopes, such 
as two neutron separation energies Son, radius, density distribution, single particle energy 
level, contribution from continua, spin-orbit splitting, potential diffuseness. In Sec. IV the 
prospect on the theoretical progress and experimental expectation of the giant halo nuclei 


are reviewed. Finally, Sec. V summarizes our main results. 


II. RELATIVISTIC CONTINUUM HARTREE BOGOLIUBOV THEORY 


The RCHB theory, which is the extension of the RMF theory with the Bogoliubov 
transformation in the coordinate representation, is suggested in Ref. [9], and its detailed 
formalism and numerical solution for a particular case can be found in Ref. [10] and the 
references therein. The basic ansatz of the RMF theory is a Lagrangian density whereby 
nucleons are described as Dirac particles which interact via the exchange of various mesons 
and the photon. The mesons are the scalar sigma (a), vector omega (w) and iso-vector vector 
rho (p). The rho (Ø) meson provides the necessary isospin asymmetry. The scalar sigma 
meson moves in a self-interacting field having cubic and quadratic terms with strengths go 


and g3, respectively. The Lagrangian then consists of the free baryon and meson parts and 


L the interaction part with minimal coupling, together with the nucleon mass M, and Mo, go, 


N Mu, Ju; Mp, Jp the masses and coupling constants of the respective mesons and: 


L = bi — M) + 48,000 — U(a) — 40,00” 
Hm ww” — hah + ymp,” — $F F” (1) 


—goow — guid — gpbpry — ep Ay. 


© The field tensors for the vector mesons are given as: 


QY = ofw” — O wt, 
Re = OM pr — o — gr x pr), (2) 
Fe = on dv — pv Au. 


~ For a realistic description of nuclear properties, a nonlinear self-coupling for the scalar mesons 


turns out crucial [17]: 
Ulo) = Imo? + 2o + Bot (3) 
(is 3 4 
The classical variation principle gives the following equations of motion : 


la-p+W(r) + BM + Vs(T)) Y: = eyi (4) 


for the nucleon spinors and 


for the mesons, where 


> 


= gb + g7 + zel- n), 
Vs(F) = goat?) 


5 
ai 
| 


(6) 
are the vector and scalar potentials respectively and the source terms for the mesons are 


Ps = 23 wii 
eR) = CA didi 
N iH (7) = Ya piy Thi 


. - „l-r: 
BA = Ch ber 


WPi, 


) where the summations are over the valence nucleons only. It should be noted that as 
N usual, the present approach neglects the contribution of negative energy states, i.e., no-sea 


) approximation, which means that the vacuum is not polarized. The coupled equations Eq. 


F 
as 


4) and Eq. (5) are nonlinear quantum field equations, and their exact solutions are very 
ø complicated. Thus the mean field approximation is generally used: i.e., the meson field 


5 operators in Eq. (4) are replaced by their expectation values. Then the nucleons move 


II 


© independently in the classical meson fields. The coupled equations are self-consistently 
~ solved by iteration. 

For spherical nuclei, i.e., the systems with rotational symmetry, the potential of the 

nucleon and the sources of meson fields depend only on the radial coordinate r. The spinor 

. f . . ; 1 

is characterized by the angular momentum quantum numbers l, j, m, the isospin t = +5 

for neutron and proton respectively, and the other quantum number 7. The Dirac spinor has 


the form: 


where Yj,,,(9, 6) are the spinor spherical harmonics and Gi (r) and FË (r) are the remaining 


radial wave function for upper and lower components. They are normalized according to 
f GEOP + FP@)P) =1. (9) 


The radial equation of spinor Eq. (4) can be reduced as : 


CEF (r) + (M + Va(r) + W(r) G(r) io 
aki (r) = (452+ 2)GP(r) (M + Valr) — Wr) F t). 


where 


= —(j + 1/2) for j =l + 1/2 
+(j + 1/2) for j =l — 1/2. 


N The meson field equations become simply radical Laplace equations of the form: 


(iE m) o= s0) (11) 


mg are the meson masses for ¢ = 0,w,p and for photon ( me = 0 ). The source terms are: 


— gops — 9207(r) — g30%(r) for the ø field 


\ Ju Pv for the w field 
tle sg(7) = (i2) 
i IpP3(r) for the p field 
a ep-(r) for the Coulomb field, 
Anrp.(r) = Dei (IGi(r))? — Flr) 
Anr?p(r) = CEi(IGi(r)? + Flr) (13) 
Amr ps(r) = Lea(lGp(r))? + FI — Sela Gar)? + |Fa(r)|?) 
4rr?pelr) = Lp-a(lGp(r)|? + |F). 
The Laplace equation can in principle be solved by the Green function: 
(r) = f rdr'Gelr,r')s4(r"), (14) 
where for massive fields 
1 1 ' , 
ipe oat me|r—r'| me|r+r’| 1 
Ganr) = Fe emia) (15) 


© It is obtained from one-meson exchange interaction Vp, pg (TT; rar 


and for Coulomb field 


p l/r forr > r 
G(r, r’) = (16) 
l/r’ forr < r. 

Eqs. (10) and (11) could be solved self-consistently in the usual RMF approximation. For 
RMF, however, as the classical meson fields are used, the equations of motion for nucleons 
derived from Eq. (1) do not contain pairing interaction. In order to have pairing interaction, 
one has to quantize the meson fields which leads to a Hamiltonian with two-body interaction. 
Following the standard procedure of Bogoliubov transformation, a Dirac Hartree-Bogoliubov 


equation could be derived and then a unified description of the mean field and pairing 


correlation in nuclei could be achieved. For the detail, see Ref. [10] and the references 


T= therein. The RHB equations are as following: 


where 


C l is the Dirac Hamiltonian and the Fock term has been neglected as is usually done in RMF. 


= The pairing potential is : 


Ane (T, P) = far, J dri > View ikl P Tri) Ray (Ti 7). (19) 
kek! 


r3) in the pp-channel and 


the pairing tensor k = V*UT 
Kaw (FP) =< Jazar] >= OE (F) bE (P) (20) 
The nuclear density is as following: 
pli. -SE a (21) 


As in Ref. [10], Vw gy used for the pairing potential in Eq. (19) is either the density- 
dependent two-body force of zero range with the interaction strength Vo and the nuclear 


matter density po: 


V (rı, T2) = Voð (rı = 2); [1 = T102] h = a) (22) 


or Gogny-type finite range force with the parameter u;, W;, Bi, H; and M; (i = 1,2) as the 

finite range part of the Gogny force [18]: 
Virar) = So eri ra)/m (W, + BP? — HiP" — M,P7P’) (23) 

i=1,2 
A Lagrange multiplier À is introduced to fix the particle number for the neutron and proton 
as N = Trp, and Z = Trpp. 

In order to describe both the continuum and the bound states self-consistently, the 
RHB theory must be solved in coordinate representation, i.e., the. Relativistic Continuum 
Hartree-Bogoliubov ( RCHB ) theory [10]. It is then applicable to both exotic nuclei and 
> normal nuclei. In Eq. (17), the spectrum of the system is unbound from above and from 
5 below the Fermi surface, and the eigenstates occur in pairs of opposite energies. When 


ç ) spherical symmetry is imposed on the solution of the RCHB equations, the wave function 


can be conveniently written as 


¡ev t) ev) 
3 =] päi oe Vil, pxl) Y PO G3 Vin(O,0)xe(t). 24 


The above equation Eq. (17) depends only on the radial coordinates and can be expressed 


¢ - as the following integro-differential equation: 


eu) K 


+ ŽGu(r) — (E +A — Wir) + Vale) Fulr) +r f r'dr'Alr, r’) Fy(r') = 0 
aio — 2 Fy(r) + (E + À — Vv(r) = Vs(r))Gu(r) +r | r'dr'Alr,r')Gv(’) = 0 
ddv) + ŽGy(r) + (B= A + Wr) = Vs(r))Fy(r) +r | r'dr'Alr, r) Fu(r') = 0 a 

d if oa ~Fy(r) — (E — A + W(r) + V(r) )Gv(r) +r I r'dr'A(r,r")Gu(r^) = 0, 


where the nucleon mass is included in the scalar potential Vs(r). Eq. (25), in the case of 
6-force given in Eq. (22), is reduced to normal coupled differential equations and can be 
solved with shooting method by Runge-Kutta algorithms. For the case of Gogny force, the 
coupled integro-differential equations are discretized in the space and solved by the finite 


element methods, see Ref. [10]. Instead of solving Eqs. (10) and (11) self-consistently for 


9 


the RMF case, now we have to solve Eqs. (25) and (11) self-consistently for the RCHB case. 
As the calculation for Gogny force is very time-consuming, we use them for one nucleus and 


fix the interaction strength of the 6-force given in Eq. (22). 


II. RESULTS AND DISCUSSIONS 
A. Binding energies and two neutron separation energies 


The numerical techniques of the RCHB theory can be found in Ref. [10] and the references 
therein. In the present calculations, we follow the procedures in Ref. [10] and solve the RCHB 
equations in a box with the size R = 20 fm and a step size of 0.1 fm. The parameter set 


NL-SH [19] is used, which aims at describing both the stable and exotic nuclei. The use of 


Š the TM1 parameter set provides similar results and we show here only those of the NL-SH 


N parameter [20]. The density dependent ð-force in the pairing channel with po = 0.152 fm~? 
C : is used and its strength Vo = 650 MeV - fm~? is fixed by the Gogny force as in Ref. [10]. 
- The contribution from continua is restricted within a cut-off energy Ea, ~ 120MeV. 

In this work the RCHB calculation is restricted to the spherical shape, which is a good 
c 2 approximation for most proton magic nuclei with Z = 8, 20, 28, 50, 82. The RCHB code 


I was also carried out for the Zr isotopes, with the considerations of the sub-magic number 


Se Z = 40 and the fact that most of the investigations with non-relativistic codes and with 
= relativistic RMF+BCS codes show that the nuclei above !**Zr are spherical [11]. 

: We calculated ground state properties of all the even-even O, Ca, Ni, Zr, Sn, Pb isotopes 
ranging from the proton drip-line to the neutron drip-line with the RCHB code. We here 
present the binding energy F» calculated from the RCHB method for these isotope chains 
and the corresponding data available [21] in the last two columns in table I-IV. We do not 
include the tables for the Sn and Pb isotopes here to save the space. The difference (AE) 
between the experimental and calculated binding energies is less than 3.5 MeV for most 
nuclei, which is less than 1% of the experimental values. The large difference (AE, ~ 10 
MeV) for some Zr isotopes are found, which are mainly due to deformation. The two neutron 


separation energy Sən defined as 
Son = lZ, A+ 2) — lZ, A) (26) 


10 


is quite a sensitive quantity to test a microscopic theory. The two neutron separation energy 
becomes negative when the nucleus becomes unstable against the two neutron emission. 
Hence, the drip line nucleus for a corresponding isotope chain is the one before the nucleus 
with the negative Sən. In Fig. 1 both the theoretical and the available experimental Sn are 
presented as a function of neutron number N for the O, Ca, Ni, Zr, Sn, Pb isotope chains, 


respectively. The good coincidence between experiment and calculation is clearly seen. 
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FIG. 1. The two neutron separation energies Sən for the proton magic isotopes are plotted as a 
function of the neutron number. The figures from left to right and from up to down are those for 
O, Ca, Ni, Zr, Sn and Pb, respectively. The experimental values are denoted by the solid symbols 


and the calculated ones by the open symbols. 


TABLE I. The ground state properties of even O isotopes calculated with the RCHB theory 
with the NL-SH parameter. Listed are neutron(r,,), proton(r,), matter(rm), and charge(re) root 
mean square radii as well as binding energies Æp. The corresponding experimental values re and 


E, are also included when available. 


Nucleus Radius (fm) Ey (MeV) 


11 


40 
120 
1409 
160 
180 
200 
20 
240 
260 
280 
300 


exp. 
c Tn Tp 
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Timi Te Exp. RCHB 
-58.549 -62.169 
-98.733 -99.641 
-127.619 -128.347 
-139.807 -142.081 
-151.370 -153.397 
-162.030 -163.124 
-168.500 -170.169 

-174.126 
-177.040 
-174.453 


Nucleus 
ACa 
34Ca 
36Ca 
38Ca, 
40 Ca 
42Ca 
Ca 
46Ca 
48 Ca, 
50Ca 
52Ca 
54Ca 
56Ca 
58Ca 
60Ca, 
62Ca 
64Ca 
66 Ca, 
68 Ca, 
Ca 
72g 


T4Ca 


TABLE II. The same as Tab. I, but for even Ca isotopes. 


Radius (fm) 


Tp 
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E, (MeV) 

Exp. RCHB 
-245.600 -244.911 
-281.360 -278.557 
-313.122 -310.051 
-342.052 -340.006 
-361.895 -360.995 
-380.960 -380.086 
-398.769 -397.937 
-415.991 -414.910 
-427.491 -424.443 
-436.600 -433.223 
-443.800 -441.354 
-449.600 -448.998 

-456.304 
-463.260 
-465.323 
-466.562 
-467.320 
-467.752 
-467.957 
-467.993 
-467.645 


Nucleus 


ANi 


3.776 
3.813 
3.842 
3.860 


TABLE III. The same as Tab. I, but for even Ni isotopes. 


3.344 
3.414 
3.476 
3.932 
3.582 
3.668 
3.745 
3.817 
3.886 
3.954 
4.025 
4.075 
4.122 
4.168 
4.211 
4.252 
4.357 
4.457 
4.552 
4.639 
4.717 
4.782 
4.835 
4.882 
4.925 
4.967 


Radius (fm) 


Tp 
3.669 
3.651 
3.641 
3.634 
3.630 
3.658 
3.686 
3.712 
3.736 
3.758 
3.778 
3.794 
3.809 
3.826 
3.842 
3.858 
3.873 
3.888 
3.902 
3.917 
3.934 
3.954 
3.979 
4.006 
4.034 
4.064 
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3.537 
3.549 
3.566 
3.585 
3.606 
3.663 
3.718 
3.770 
3.821 
3.872 
3.925 
3.965 
4.003 
4.042 
4.079 
4.115 
4.194 
4.271 
4.346 
4.417 
4.482 
4.540 
4.592 
4.638 
4.683 
4.727 


3.755 
3.738 
3.727 
3.721 
3.717 
3.745 
3.772 
3.797 
3.821 
3.842 
3.862 
3.877 
3.893 
3.908 
3.924 
3.940 
3.955 
3.969 
3.983 
3.998 
4.014 
4.035 
4.058 
4.085 
4.113 
4.142 


E, (MeV) 

Exp. RCHB 
-348.814 

-385.500 -384.925 
-420.460 -419.044 
-453.150 -451.786 
-483.988 -483.498 
-506.454 -503.244 
-526.842 -522.146 
-545.259 -540.336 
-561.755 -557.802 
-576.830 -574.370 
-590.430 -589.508 
-602.600 -602.094 
-613.900 -613.412 
-623.900 -623.982 
-633.100 -634.009 
-641.400 -643.615 
-646.356 

-648.979 

-651.488 

-653.855 

-656.036 

-658.026 

-659.852 

-661.570 

-663.220 

-664.767 


100.7; 
Ni -663.123 
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Nucleus 
AZr 
T4Zr 
T6Zr 
78 Zr 
80 Zr 
82 Zr 
ca A 
36 Zr 
88 Zr 
Zr 
9272r 
%4Zr 
36 Zr 
98 Zr 
100 Zr 
10277 
1047, 
106 Zr 
108 Zr 
110 Zr 
1127, 
1147, 
116 Zr 
118 Zr 
120 Zr 
12277 


1247r 


126 Zr 


4.270 
4.305 
4.330 
4.349 


TABLE IV. The same as Tab. I, but for even Zr isotopes. 


3.932 
3.984 
4.035 
4.088 
4.133 
4.174 
4.214 
4.251 
4.286 
4.346 
4.403 
4.458 
4.511 
4.563 
4.613 
4.664 
4.716 
4.769 
4.817 
4.856 
4.889 
4.920 
4.951 
4.981 
5.010 
5.103 
5.192 


Radius (fm) 


Tp 
4.149 
4.150 
4.155 
4.163 
4.166 
4.168 
4.171 
4.176 
4.181 
4.204 
4.226 
4.248 
4.269 
4.288 
4.306 
4.323 
4.339 
4.355 
4.370 
4.385 
4.401 
4.417 
4.433 
4.449 
4.466 
4.474 
4.483 
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4.051 
4.072 
4.097 
4.126 
4.149 
4.171 
4.194 
4.217 
4.240 
4.285 
4.329 
4.372 
4.414 
4.455 
4.496 
4.536 
4.578 
4.620 
4.660 
4.693 
4.723 
4.753 
4.782 
4.810 
4.838 
4.909 
4.977 


4.225 
4.226 
4.231 
4.239 
4.242 
4.244 
4.247 
4.252 
4.257 
4.280 
4.302 
4.323 
4.343 
4.362 
4.380 
4.397 
4.413 
4.428 
4.443 
4.458 
4.473 
4.488 
4.504 
4.520 
4.537 
4.545 
4.553 


Ey (MeV) 

Exp. RCHB 
-568.120 

-600.450 

-631.267 

-669.800 -659.843 
-694.700 -686.399 
-718.190 -711.679 
-740.640 -736.106 
-762.606 -759.890 
-783.893 -783.173 
-799.722 -795.603 
-814.677 -807.665 
-828.994 -819.421 
-840.972 -830.903 
-852.440 -842.122 
-863.720 -853.053 
-874.500 -863.613 
-884.000 -873.647 
-892.300 -883.105 
-891.983 

-900.136 

-907.646 

-914.714 

-921.459 

-927.963 

-934.285 

-934.256 

-934.271 


128 Zr 
130 Zr 
132 Zr 
34Zr 


136 Zr 


387r -934.540 


5.275 
5.354 
5.429 
5.501 
5.563 


4.491 
4.499 
4.508 
4.516 
4.527 
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5.043 
5.106 
5.167 
5.226 
5.280 


4.562 
4.570 
4.578 
4.587 
4.597 


-934.329 
-934.419 
-934.526 
-934.618 
-934.640 


We See from Fig. 1 and Tab. I-IV the proton magic even mass nuclei at the neutron 
drip-line are predicted as 730, "Ca, SNI, °Zr, '4Sn, and ?°°Pb, respectively. Out of these 
neutron drip line nuclei, we have experimental information only for the O isotope. The 
experimental efforts were made to investigate 7°O [22,23] and *8O [24] whether they are 
bound. These nuclei were found unstable and the heaviest O isotope was concluded to be 
240. Hence, the present calculation is not successful in reproducing the O neutron drip line. 
Numerous theoretical studies of binding energies of the neutron-rich O nuclei have been 
made [25], but these calculations failed to reproduce that ™*O is the drip line nucleus as the 
case of RCHB with NL-SH presented here. It is known that in the relativistic mean field 
theory, one parameter set is used to describe all nuclei in the nuclear chart, which is quite 
a challenge. In fact, it is well known that the mean field theory is difficult for light nuclei. 

ka Therefore we have NL1, NL-SH, NL3, and TM1 for the heavy system and NL2 and TM2 
N = for the light system. The heavy oxygen isotopes lie just at the border line between the light 
© and heavy groups. Therefore it is not surprising that the NL-SH, NL1, NL2, NL3, TM1, 


se, TM2 and the non-relativistic counterparts mostly do not predict correctly the neutron drip 
<= line which is already known experimentally to be at 740. 


One may suspect that the drip-line nucleus for Ca is predicted to be Ca instead of 


` ~ Ca. Most calculations predict the Ca isotopes with N > 50 would be unbound, as the 
= more neutrons should fill the continuum region above the N = 50 shell. For example, in 


© Ref. [13] and Ref. [14], the drip-line nucleus is predicted as Ca with the HFB and Skyrme 


‘= HF method, respectively. 

© The reason for the unexpected bound nuclei ‘Ca is caused by the halo effect, which also 
results in the disappearance of the normal N = 50 magic number in the RCHB calculation. 
Seen along the Sən to N curve in Fig. 1, strong kinks can be clearly seen at N = 8, 
N = 20, 28,40, N = 28,50, N = 50,82 and N = 126 for O, Ca, Ni, Zr, Sn, Pb isotope 
chains, respectively. All these numbers correspond to the neutron magic numbers, which 
come from the large gaps of single particle energy levels. However, no kink appears at the 
N = 50 magic number for Ca isotopes in Fig. 1. In stable nuclei, the N = 50 magic 
number is given by the big energy gap between the 1g9/2 and 3s, 2d levels. However in 


the vicinity of the drip-line region of Ca isotopes, the single particle energy of 3s,/2 would 
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decrease with the halo effect and its zero centrifugal potential barrier, the large gap between 
1g9/2 and 3s1/2 disappears. Therefore, the nucleus ‘Ca is no more a double-magic nucleus 
and furthermore 'Ca is also bound. This disappearance of the N = 50 magic number at 
the neutron drip line is due to the halo property of the neutron density with the spherical 
shape being kept, different from the disappearance of the N = 20 magic number due to 
deformation (the unbound 7%?8O mentioned above). Most recently a new magic number 
N = 16 has been discovered in neutron drip-line light-nuclei region which is also considered 
as owing to the halo effect [4]. It can be expected also that near drip line, more disappearance 
of the traditional magic numbers and regeneration of new magic number would be found 
with the same mechanism. 
Another remarkable phenomenon in Fig. 1 is that the Sən values for exotic Ca iso- 
R j topes near neutron drip line are very close to zero in the large mass region, i.e., Som 7% 


~N 2.06, 1.24, 0.76, 0.43, 0.21, 0.04 MeV for A = 62, 64, 66, 68, 70, 72 isotopes, respectively. 


) If one regards °°Ca as a core, then the several valence neutrons occupying levels above the 
oð N = 40 sub-shell for these A > 60 exotic nuclei are all weakly bound and can be scattered 
= easily to the continuum levels due to the pairing interaction, especially for ĉ6-72Ca. This 
p case is very similar to Sọ» in Zr isotopes with N > 82 (see Fig. 1) [11] . Note that for 
Sn isotopes in the vicinity of the drip line [26], however, Sən decrease rather fast with the 
= mass number and for Ni isotopes [10], S2n are quite large (~ 2MeV) and undergo a sudden 
© drop to negative value at the drip-line nucleus. Such behavior of Son in Ca isotopes with 
"= N > 40 gives us a hint that there exist giant halos in these nuclei, just as what happens in 


© Zr isotopes with N > 82 [11]. 


B. Nuclear root mean square radii 


We discuss here the nuclear radii, which are the important basic physical quantity to 
describe atomic nuclei as well as the nuclear binding energies. In the mean field theory, the 
root mean square (rms) proton, neutron, mass radii (rp, Tn, Tm) can be directly deduced 
from their density distributions p: 


2 1/22 (L Perdar 1/2 (27) 


—_ T opdr 
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and 


2 
=e r2, s (Le + pn)r°dr 1/2, (28) 
J (Pp + fn)dt 
The theoretical rms charge radii re are connected with the proton radii rẹ, as: 
r= r + 0.64 fm?, (29) 


when the small isospin dependence is ignored for exotic nuclei [27]. 
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FIG. 2. The root mean square charge radii for the proton magic even mass nuclei are plotted as 


" a function of the neutron number. The dotted curves correspond to the radii calculated with the 


A*/3 formula. 
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FIG. 3. The root mean square neutron radii for the proton magic even mass nuclei are plotted 
as a function of the neutron number. The solid curve corresponds to the radii calculated with the 


N*/3 formula. 
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FIG. 4. The root mean square radii for the proton, neutron, charge and matter distributions are 


w 
O 


plotted as a function of the neutron number for the Ca isotope. The available experimental values 


are denoted by solid symbols. 


— It is interesting to investigate the neutron number dependence of the proton radii (or 


[© charge radii) and the neutron radii for these proton-magic isotope chains. In Fig. 2, we 


| show the rms charge radii re obtained from the RCHB theory (open symbols) and the data 


) available (solid symbols) for the even-even Ca, Ni, Zr, Sn and Pb isotopes. The values for 


the O isotopes are not plotted here considering the systematic deviation of the mean field 
|> theory for light isotopes as discussed above. As it could be seen, the RCHB calculations 
© reproduce the data very well (within 1.5%). For a given isotopic chain, an approximate 


5 linear N dependence of the calculated rms charge radii re is clearly seen. However, the 


>< variation of re for a given isotopic chain deviates from the general accepted simple A'/? law 
7 (denoted by dashed lines in Fig. 2), which shows a strong isospin dependence of nuclear 
= a charge radii is necessary for nuclei with extreme N/Z ratio. Recently, we have investigated 
i these charge radii as well as the large amount of experimental data and extract a formula 
with the Z'/? dependence and isospin dependence to better describe these charge radii all 
over the nuclear isotope chart [28]. 

We also plot in Fig. 3 the neutron radii r,, from the RCHB calculation for the even-even 
nuclei of the O, Ca, Ni, Zr, Sn and Pb isotope chains. The prediction r, curve using the 
simple experiential equation rn = ro- N‘/? with ro = 1.139 normalizing to 2°°Pb is also 
represented in the figure. Except for the O and Ca isotopes, it is very interesting to see that 
the simple formula for r, agrees with the calculated neutron radii except for two anomalies. 


One appears in the Ca chain above N = 40, and the other in Zr chain above N = 82. The 
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increase of exotic Ni and Sn nuclei is not as rapid as that in Ca and Zr chains. The regions 
of the abnormal increases of the neutron radii are just the same as those for Sn. Both 
behaviors are connected with the formation of giant halo. As Calcium is a light element, the 
giant halos in exotic Ca isotopes would be much more easily accessed experimentally than 
in the Zr chain. 
In Tab. II we present the calculated neutron, proton, matter, and charge radii for all even 
Ca isotopes, as well as several available charge radii data. Those radii against mass number A 
are plotted in Fig. 4. For nuclear charge radii re, the well-known parabolic behavior along 
40-48Ca is not reproduced in our calculation. It is mainly due to the improper spherical 
supposition to the real deformed *:**“°Ca nuclei. It can be clearly seen that the radii rp, 
Tn, Tm and re all increase with the mass number A. The mass radii rm as well as r, increase 
a much faster than rp and re. Besides the normal increase of r, and rm with A, a slightly 


~ up-bend tendency occurs in the neutron drip line region from A = 62 up to A = 72. 


C. Density distribution and halos 


We shall concentrate from here the giant halo properties for the Ca isotopes as mainly 
| due to the easier access by the present experimental techniques. For exotic nuclei, physics 
connected with the low density region in the tails of the neutron and proton distributions 


` have attracted a lot of attention in nuclear physics as well as in other fields such as astro- 


if 


physics. It is therefore of great importance to look into the matter distribution and see 
‘© how the densities change with the proton to neutron ratio in these nuclei. As the density 
here is obtained from a fully microscopic and parameter free model which is well supported 
by the experimental binding energies, we now proceed to examine the density distributions 
of the whole chain of Calcium isotopes in this section and study the relation between the 


development of halo and shell effects within the model. 
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r (fm) 


FIG. 5. (a) The neutron, (b) proton, and (c) matter density distributions for the Ca isotopes. 


The arrows represent the density variation tendency with the neutron number JN. 


The density distributions for some Ca isotopes are given for the neutron, proton and 
» matter densities in Fig. 5a-c, respectively. In these figures, the arrows represent the density 
variation tendency with the mass number A. From Fig. 5c, the matter densities extend 
outwards with A at the nuclear surface (3 < r < 6 fm), while the variation is small in the 
center (r < 3 fm). For the proton density in Fig. 5b, the density for r < 3 fm decreases 
dramatically, and at the surface (r ~ 4fm) extends outwards monotonously with A, although 
the proton number is fixed to 20. It is in accordance with the increase of the proton radii 
Tp and the charge radii re as shown in Fig. 4. The increase of the radius is just due to the 
nuclear saturation property, as the total mass number A of the nuclei is increasing. The 
neutron densities, which are represented in Fig. 5a, extend greatly outwards with A on the 


surface and neutron skin is developing. In the center, most neutron densities also increase 
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with A, however the increase is not as dramatic as the decrease of the proton density. As a 
result, the matter densities vary (mainly decrease) slightly in the interior region. 

The halo phenomena are always related with the wide extension of the nuclear density 
distributions in space. In order to exhibit clearly the large spacial extension, the density 
distributions are also shown on a logarithmic scale in the corresponding figures of Fig. 5a, 
b, c as the inset for neutron, proton and matter density distribution. From the inset figures 
of Fig. 5a, the proton density decreases once more with A in the rather exterior region of 
r > 8 fm, meanwhile heavier Ca isotopes extend toward outside in the region of 4 < r < 6 
fm and the proton density decreases dramatically with A in the center r < 3 fm. So there 
exist two switching regions for the proton density: one lies in the nuclear surface (r ~ 4 fm) 


due to the nuclear saturation property, and the other appears over 6 — 8 fm region, which 


' are related with the neutron-deficiency. 

l From the neutron density represented in the inset of Fig. 5a, it is clearly seen that 
' the neutron density distributions for Ca isotopes with N > 40 extend much more widely 
» than those for lighter ones. The tail becomes larger and larger with the neutron number N 
— increased from 42 to 52, which is quite different from the cases in Ni and Sn isotopes [10,26]. 
, In the latter two isotope chains, the tails of the neutron densities reach saturation beyond 


; a specific nucleus, i.e., A = 90 for Ni and A = 160 for Sn. However, no saturation point is 


; seen for the Ca isotopes as shown in Fig. 5a. 

For the tails of the matter density distributions in the inset of Fig. 5c, they are mainly 

determined by the proton densities for proton-rich nuclei while by the neutron densities for 

“ neutron-rich nuclei. The nucleus with the smallest tail for matter distribution is shown to 
be “Ca, which is nearest to the 6-stability line, while not the lightest neutron-deficiency 
nucleus °4Ca. For neutron rich nuclei with N from 42 to 52, the tail behaviors of rm are 
consistent with that of r,, which are strong evidences that giant neutron halos exist in exotic 


Ca nuclei (A > 60). 
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FIG. 6. The radii r corresponding to different given proton and neutron densities 
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Pp(T)(Pr(r)) = 107?,1074, and 10~° fm~3),respectively, as a function with the mass number in 


Ca isotope chain. 


To give a simple quantitative idea on the density distributions and halos, in Fig. 6, the 


~ radii at which the proton and neutron densities, p,(r) and pp(r) = 107?, 1074, 1078 fm~®, 


© respectively, as a function of the neutron number for the Ca isotope chain are shown. As the 


; central density of nuclear matter is about pp ~0.16 fm~? and the proton and the neutron 
— density is about a half of pp, the value r for ppin) = 107? fm? corresponds to the densities 
™ decrease to 10% of the central density and the surface of the nuclei. The value r for ppn) = 
EN 1076 fm™? corresponds to the tails of density distribution. From this figure, for *7Ca, we see 
= that the proton and neutron densities have the same values at the same r, i.e., the neutron 
€ 5 density distribution is very similar to that of proton for the 6 stable nuclei. For neutron 
= rich nuclei, if pp = pp, the neutron radius r is much bigger. For neutron deficient nuclei, 
© the opposite are seen. For pp = 10~° fm~3, a increasing tendency for r curve can be clearly 

seen, particularly, the curve shows a strong kink at N = 62. The kink is also a important 


signal for the giant neutron halos, as well as the Son and neutron radii rẹ. 


D. Single particle levels in canonical basis and contribution of continuum 


The above results can be understood more clearly by considering the microscopic struc- 
ture of the underlying wave functions and the single-particle energies in the canonical basis 
[10]. In Fig. 7, the neutron single-particle levels in the canonical basis are shown for the 


even Ca isotopes from the mass number A = 34 to A = 72. The shell closure (N = 20 and 
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28) and sub-shell closure (N = 40) are clearly seen as big gaps between levels. A dotted 
line in the figure represents the neutron chemical potential An, which jumps three times at 
the magic or submagic neutron number on its way to almost zero at ‘Ca. These jumps 
correspond to the shell (or sub-shell) closure just as the appearance in the Sə» case. The 
An comes close to zero for nuclei near the neutron drip line, i.e., An œ% —0.92 MeV for °Ca, 
-0.64 MeV for Ca, -0.44 MeV for Ca, -0.28 MeV for Ca, -0.16 MeV for Ca, and -0.07 
MeV for “Ca. Meanwhile, no jump at N = 50 in the \,, curve is in agreement with the 
disappearance of this traditional magic number mentioned above in the Sn case. From Fig. 
7, the big gap between the lgg/2 orbit and the s — d shell has disappeared in the neutron 
drip line region due to the lowering of 31/2 and 2d5/2 orbits. For example, this gap for Ca 


is only 1.02 MeV. 


10 m~ SF 
ae 19o/0 2d;z 


Single Particle Energy [MeV] 


FIG. 7. The single-particle neutron energies in the canonical basis with the mass number A for 


? the Ca isotopes. The Fermi surfaces are shown by the dot-dashed line. 


As neutron chemical potential A, approaches zero, the Ca isotopes with A > 60 are 
all weakly bound nuclei. This means that the additional neutrons will occupy the weakly 
bound states, which are very close to the continuum region. These neutrons supply very 
small binding energies, and result in nearly vanishing two neutron separation energies Son. 
Furthermore, the pairing interaction scatters the neutron pairs from such weakly bound 
states to continua as the Fermi level is close to zero. The single particle occupation of 
continua then affects the nuclear properties. For Ca isotopes with A > 60, the added 
neutrons will occupy the weakly bound states and continua: 1499/2, 351/2, 2d5/2, 2d3/2, etc. 


Such orbits play importance roles as discussed below. 


26 


-20 -15 -10 -5 -15 -10 -5 0 
E [MeV] 


FIG. 8. The occupation probabilities in the canonical basis for various even Ca isotopes as a 


function of the single particle energy. The chemical potential is indicated with a vertical line. The 


number Np of neutron in the halo are also shown. 


In Fig. 8, we present the occupation probabilities v? of neutron levels near Fermi surface 
(i.e., —20 < E < 10 MeV) in the canonical basis for several neutron-rich even Ca isotopes, 
> i.e., 58-72Ca. The neutron chemical potential \,, is indicated by a vertical line. For nuclei 
with mass number A < 60, the chemical potential is quite large (e.g., -13.1 MeV for Ca, 
-6.74 MeV for *8Ca, and -3.69 MeV for *°Ca), the occupation probabilities for the continua 
are nearly zero. As the neutron number goes beyond the subshell N = 40, the Fermi 
surface approaches to zero and the occupation of the continuum becomes more and more 
important. Summing up the occupation probabilities v? for the states with E > 0, one can 
get the contribution of the continua n, [11]. They are approximately 2.2, 0.6, 1.1, 1.7, 1.9 
and 2.7 for Ca, ®Ca, Ca, Ca, Ca, and Ca, respectively. In ref. [11], from two to 
roughly six neutrons in total scattering in the continuum is predicted in Zr isotopes with 
N > 82, while the neutron number scattering to continuum in Ca isotopes is reduced to 


about 1 ~ 3. 
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FIG. 9. The single-particle levels in the canonical basis for the neutron in ®°Ca. The neutron 


© potentials V(r)+ S(r) is denoted by the solid line and the Fermi surface is shown by a dashed line. 


j = The occupied probabilities for the single-particle levels are shown by their line-length. Also the root 
) mean square radius rj in fm for every level is listed in the parenthesis behind the corresponding 


level. 


As a typical example, we investigate in detail the single particle levels and contributions 


5 < of continua with exotic nuclei ©°Ca. Fig. 9 shows the neutron single-particle levels in the 


! canonical basis for Ca. The length of each level is proportional to its occupied probability 
v’. The Fermi surface for neutron (A, œ~ —0.435 MeV) is represented with a dashed line, 
and the nuclear mean field potential V(r) + S(r) is denoted with the solid curve. We here 
define the root mean square radius r,,; for one single-particle level denoted by nlj as: 
J Pnig? 2dr 


; 30 
Í PajdT ( ) 


Talj = 


where prj presents the probability density of each level and the corresponding wave function 
Wnij- This root mean square radius fm; for each level is shown in the parenthesis behind the 
corresponding level label in units of fm. Also in Tab. V, we list the single particle energies 
€nij, the rms radii ruj and the occupied probabilities viij for all the neutron levels with 


energies € < 10 MeV in Ca calculated from RCHB theory. 
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Here, the state 1g9/2 is weakly bound with the energy -0.48 MeV, while the 351/2, 2d5/2, 
2d3/2, and 1g7/z levels are in the continuum with energies 0.64, 1.41, 2.85, and 5.67 MeV, 
respectively. We note that due to the absence of a centrifugal barrier, the orbit 3s1/2 lies 
below states 2d3/2 and 2ds/2. The occupied probabilities v? of these states are: 0.514 for 
1g9/2, 0.089 for 31/2, 0.071 for 2ds5/2, 0.027 for 2d3/2, and 0.017 for 1g7;2. We get about 
0.85 neutrons in these continuum states, which are about 4/5 of the total contribution of 
the continuum np = 1.1 for Ca (see Fig. 8). The 31/2 state has a rms radius 7.24 fm in 
comparison with the rms radii of the neighbor states (~ 5 fm) and the total neutron rms 
radii (4.314 fm) owing to its zero centrifugal barrier. Therefore the nucleon occupying on 
the state 351/2 will contribute to the nuclear rms radius considerably. 

The relative contributions pyij;/p, of different single-particle orbits to the full neutron 
density as a function of radical distance r are represented in Fig. 10 for the typical nuclei 
66Ca. For comparison we also present the total neutron density p, with the shaded area in 
) arbitrary units. In the interior of nuclei (r < 4 fm), the contribution to neutron distribution 
» mainly comes from the low-lying energy levels such as 151/2, 1p3/2, lpi/2, 1d5/2, 281/2, 1d3/2. 
In the nuclear surface(r ~ 5 fm) and the beginner part of the tail (6 < r < 10 fm), the 
ə weakly bound f-p shell states and very weakly 1go/2 orbit play dominant role for neutron 


\ distribution. Then their contributions gradually decrease with the radius, while the contri- 


“== butions from the continuum (e.g. 381/2, 2d5/2) become more and more important. For r > 15 


5 fm, the state 3s; /2 is dominant (about 60 percent) owing to its zero centrifugal barrier, and 
the other continuum states 2d5/2, 2d3/2,481/2 and the very weakly bound state 1g9/2 also play 
their roles. It is well shown that the contributions of the continua are crucial to the tail, 


which is closely connected with the halo phenomena. 


TABLE V. The properties of neutron single-particle levels in ©°Ca, including single particle 


energies, occupation probabilities, and rms radii for the single-particle levels. 


€ v r 
0.64 0.089 7.24 
1.41 0.071 6.10 
2.85 0.027 6.01 
5.67 0.017 5.20 
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281/2 -20.11 1.00 4.08 481/2 7.88 0.002 8.42 


idg -20.04 1.00 3.85 3p3/2 8.11 0.002 6.51 
if -12.05 1.00 4.39 3p1/2 8.35 0.002 6.63 
2p3/2 -7.24 1.00 5.03 
1fs/2 -6.07 1.00 4.58 
212 -5.79 1.00 5.19 
199/2 -0.48 0.514 5.03 
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In Tab. VI, the energies, rms radii and occupied probabilities of the neutron single- 
particle level 31/2 for even isotopes ®°-7?Ca isotopes are presented. Its energy decreases 
with neutron number N, i.e., it lies in continuum for °°-®8Ca and becomes slightly bound 
for Ca. As both the rms radii and v? for the state 35/2 increase monotonously with 
the mass number A, it will drive the tail of neutron density distribution become larger and 
larger as shown above. Of course, with more neutrons added to Ca up to the heaviest 
bound nucleus predicted as Ca, the contributions to the tail from the other continuum 


2d5/2, 451/2 and 2d3/2 will also become more important due to their larger occupation. 
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FIG. 10. Relative contributions of the different orbits to the total neutron density as a function 


of the radius. The shaded area indicates the total neutron density in arbitrary units. 


TABLE VI. The properties of 3s1/2 for OO 2a 


A 60 62 64 66 68 70 72 


e (MeV) 3.751 1.842 1.162 0.638 0.247 -0.028 -0.238 
r (fm) 5.524 6.356 6.781 7.240 7.709 8.062 8.212 


v? ~0.0 0.015 0.039 0.089 0.201 0.404 0.628 


31 


E. Vector and scalar potentials and their isospin evolution 


One of the essential differences between a non-relativistic single particle equation and the 
relativistic Dirac equation in nuclear physics is the fact that the relativistic equation contains 
from the beginning two potentials V and S, which have different behaviors under Lorentz 
transformation [7]. We now investigate the vector and scalar potentials in the Ca isotope 
chains in detail. Fig. 11 presents the vector V and scalar potential S for proton (left panels) 
and neutron (right panels) in Ca chains as a function of the radial distance r, respectively. 
We can see from the figure that both the vector and scalar potentials have a Woods-Saxon 
shape, i.e., the similar nuclear density distribution. They nearly vanish outside the nucleus 
and they are more or less constant in the nuclear interior, namely S ~ —420 MeV as an 


attractive potential and V ~ 350 MeV as an repulsive potential. 


Proton Neutron 


V(r)/MeV 


S(r)/MeV 


r (fm) r (fm) 
FIG. 11. The vector and scalar potentials (V(r),S(r)) for the Ca isotopes. The left panel presents 


those for proton, while the right panel for neutron. 


We can clearly see that the scalar potential S(r) for proton is just the same as that for 


neutron in all cases. It is isospin independent as the scalar potential comes from ø meson 
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field, which takes the same form for either proton or neutron. On the other hand, the 
vector potential V(r) is isospin dependent: the potential for proton is somewhat different 
from neutron. The difference comes from the p meson and the Coulomb interaction. For a 
given nucleus, the vector potential for proton is a little larger than that for neutron in the 
nuclear interior, and it does not vanish outside the nucleus due to the long range Coulomb 


interaction. 


F. Spin-orbit splitting in exotic nuclei 


An advantage of the relativistic mean field theory is to get the spin-orbit coupling nat- 
urally. The spin-orbit splitting has been discussed through various themes with relativistic 
'S or non-relativistic microscopic theory [26,29,30]. Now we examine the spin-orbit splitting 


M for the whole even Ca isotopes in the relativistic continuum Hatree-Bogoliubov theory. 


, The spin-orbit splitting energy Es for the two partners (nlj = 1 —1/2,nlj = l+ 1/2) is 
4 defined as: 


Enij=1-1/2 — Enij=1+1/2 


E= 
l 2+1 


(31) 


C J In Fig. 12, the spin-orbit splitting energies E;, in Ca isotopes are shown as a function of 


= mass number A for the proton (lower panel) and neutron (upper panel) spin-orbit partners 


— 


1d3/2, 1ds/2), (1497/2, 199/2), (1p1/2, 1p3/2), (1 fs/2, Lf7/2) and (2p1/2, 23/2). For the lighter 


Il 


Ca isotopes, some doublets sit in the continuum region with large positive energies in the 
O canonical basis. It would result in some uncertainty in spin-orbit splitting. Therefore we limit 
the level energies E below 10 MeV. In Fig. 12, the splitting decreases monotonically from 
proton-rich side to the neutron rich-side for most doublets except the doublets (1p1/2, 1p3/2) 
and (2pi/2,2p3/2). The splitting of the p doublets between two closures, i.e.,“°Ca, “Ca, 
fluctuates quite a lot. Above N = 28, the splitting increases a little, then decline with the 
neutron number. For the doublets (1p1/2, 1p3/2) and (1d3/2, 1d5/2), the spin-orbit splitting 
for neutron and proton is very close to each other. For the doublets (1fs5/2,1f7/2) and 
(2p1/2, 23/2), there are quite big difference. Moreover, the splitting for neutron is usually 


smaller than that for proton. 


33 


9 © N1P si 3/2 

u—aNid 
N 1 feo, 7/2 
N1 972, 9/2 


G N2p; 2, 3/2 


3/2, 5/2 


E,, [MeV] 


Se »6-C-Ge8 -0-0 


G © PIP ip s2 
u—a Pid 


3/2, 5/2 


5/2, 7/2 


ee P2D, 10 3/2 


30 40 50 60 70 80 
A 


FIG. 12. The spin-orbit splitting energies EF), in Ca isotopes as a function of mass number A for 
the proton and neutron spin-orbit doublets (1d3/2, 1d5/2), (197/2; 199/2), (1p1/2; 13/2), (1fs/2, Lf7/2) 


and (2p1/2, 2p3/2)- 


It would be very helpful to examine the origin of the spin-orbit splitting in the Dirac 
5 equation. For the Dirac nucleon moving in a scalar and vector potentials, its equation 
of motion could be decoupled and reduced to either the upper component or the lower 
component, respectively. If it is reduced to the lower component, it will be related with 
another interesting topic — the pseudo-spin symmetry [31]. Here we pay attention to the 
spin-orbit splitting, which is related to the upper component. The Dirac equation can be 


reduced for the upper component as follows, [26]: 


ea = oe ee 
dr? E+2M -W +Vs dr dr’ 
7 we +k) 1 Kd(2M —VW + Vs) iG (r) 
r? E+2M -Ww +Vsr dr ‘ 
= —(E4+2M —W + Vs)(E —W —Vs)G2(r), (32) 


where 
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—l-1, j=14+1/2 
l j=l-1/2 


The spin-orbit splitting is provided by the corresponding spin-orbit term [26]: 


1 kd(2M - V(r) + S(r)) 


EB+2M—Vi+S@r dr 2) 


It can be seen clearly that the spin-orbit splitting is energy dependent and depends also 
on the derivative of the potential 2M — V(r) + S(r) as well as the particle distribution. 
Therefore the so-called spin-orbit potential in the RCHB theory is defined [26]: 

_ 1d(22M —V(r)+S(r)) 


t dr 


Vis 


dV (r) — S(r)) 
dr 
) Ca isotopes is given in Fig. 13. The potential difference, V(r) — S(r), for both proton and 


The derivative of the potential difference, , for proton and neutron in the 
~= neutron are almost the same, as V(r) — S(r) is a big quantity ( ~ 700 MeV ), and the 
> difference in the spin-orbit potential for proton and neutron could be neglected. Therefore 


the proton and neutron V;, are almost the same in the present model. That is the reason 


SN why the spin-orbit splitting for neutron and proton doublets (1p1/2, 1p3/2) and (1d3/2, 1d5/2) 


= is very close to each other as shown in Fig. 12. 
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FIG. 13. The derivative of the potential difference, aS: 


(right) in the Ca isotopes as a function of the radius. 
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From *“Ca to “8Ca, the amplitude of Vj, increases monotonically and from *8Ca to "Ca, 
the amplitude decreases monotonically due to surface diffuseness. Meanwhile, the maximum 
point of the potential Vs has an outwards tendency. Thus the systematic decrease of spin- 
orbit splitting is partially related with the decrease of Vis. Furthermore, the systematic 
decrease also comes from the diffuseness of the nuclear potential. In Figs. 14 and 15, we see 
that the potentials V(r) — S(r) and V(r) + S(r) extend outwards into the nuclear surface, 


which makes the diffuseness increase with the neutron number. 


Proton Neutron 


V(r)+S(r) [MeV] 


r [fm] r [fm] 
FIG. 14. The summed mean field potential, V(r) + S(r), for proton and neutron in even Ca 


isotopes. Directions of arrows in the figure show the diffuseness of nuclear potential. 
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FIG. 15. The mean field anti-nucleon potential, V(r) — S(r), for proton and neutron in even Ca 


isotopes. Directions of arrows in the figure show the diffuseness of nuclear potential. 
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FIG. 16. The spin-obit splitting Fıs as a function of a specific quantum number in 47°?:6:72Ca, 


J The left panel shows the splitting for p orbits with the quantum number n; the right panel shows 


la S the splitting for n = 1 orbits with the quantum number l. 


So far we have investigated how the proton and neutron spin-orbit splitting changes with 


™— the mass number A, i.e. the isospin dependence. The decreasing tendency is mainly caused 


€ Y by the derivative of anti-nucleon potential and the diffuseness of the mean field. For one 


= special nucleus, the potential is the same for different spin-orbit doublets, so it is necessary 


f 5 to investigate how the splitting changes with different quantum numbers, as in Fig. 16. 


= To avoid the confusion caused by showing too many doublets in the figure, we present the 


O splitting energies for Fıs of (p1/2, P3/2) partners as a function of quantum number n in the 


left panel, and those of different n = 1 doublets as a function of quantum number / in the 
right panel, respectively, for 4526272 Ca, 

In the left panel of Fig. 16, the spin-orbit splitting in doublets (1p1/2, 1p3/2) is much 
larger than that in (2p;/2, 2p3/2) cases, and that in (3p1/2, 3p3/2) for Ca is the smallest. It 


can be clearly seen in Eqn.(32) that, besides the spin-orbit potential, the spin-orbit term 
1 

E+2M —V(r)+S(r) 

E dependence. In fact, the energy E (varying from -40 to 10 MeV) is a small quantity 


compared with the potential 2M — V(r) + S(r), which is 2M (~ 1800 MeV) in surface 


also depends on the factor, However, this factor has little energy 
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or about 1000MeV in the nuclear interior. Thus the energy dependence of the spin-orbit 
splitting difference caused directly by the energy factor can be neglected. 

It has been suggested that the spin-orbit splitting difference mainly comes from the 
overlap between the density distribution and the spin-orbit potential Vis, as demonstrated 
for Sn isotopes in Ref. [26]. In Fig. 17, for “Ca, Ca, Ca and "Ca, the derivative of 
the potential V(r) — S(r), d(V(r) — S(r))/dr, and the density distribution of the p-waves 
are given in the respective lower panel, while the spin-orbit potential, Vs, multiplied by the 
density distributions for pı/2 are respectively given in the upper panel. The overlap between 
the spin-orbit potential Vi, and the particle distribution is represented by the curve in the 
corresponding upper panel. Their contribution to the spin-orbit splitting is proportional to 
_ the area surrounded by the corresponding curve and the z—axis. It is clearly seen that the 
* overlap of 1p, /2 is larger than that of 2p 1/2 for all these isotopes and that of 3p,/2 in Ca is 
the smallest. The overlap values of 2p;/2 for °*°”’*Ca are very close from each other and all 


are smaller than that for 47Ca. These observations can explain the features observed in the 


oð tight panel of Fig. 16. 
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FIG. 17. The overlap between the spin-orbit potential and density distributions for 4?Ca, °?Ca, 
— 62Ca and "Ca, respectively. In each nucleus, the spin-orbit potential Vj, multiplied by the density 
5 "4 distribution of the p;/2 orbit is given in the upper panel; the derivative of anti-nucleon potential 


- and the density distributions for pı /2 orbit in arbitrary unit is given in the lower panel. 


Seen from the right panel of Fig. 16, the spin-orbit splitting for n = 1 states has similar 
features for 4*Ca, Ca, °*Ca and “Ca. The splitting of the doublets (d/2, d3/2) is the largest 
one. With l, the splitting for n = 1 increases at first (from p to d) and then decreases (from 
d to f, g). This feature can also be understood with the overlap of the spin-orbit potential 
and the density distribution of these orbits. Here we just choose "Ca as an example (Fig. 
18), similar patterns also appear in other Ca isotopes. The largest overlap happens for 


doublets (d1/2,d3/2) in Ca is clearly seen from the upper panel of Fig. 18. 
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FIG. 18. The same as Fig. 17, but for n = 1 orbits in Ca. 


IV. THE PROSPECT OF GIANT HALOS 


Nuclear halo phenomena have been studied and investigated by many nuclear theorists 


cy and experimenters. More and more halo nuclei, including neutron halo, proton halo, and 


» halo at excited states, have been identified and reported by different methods and advanced 


“S instruments. However, only nuclei with one or two halo nucleons are found by experiment 


É until now. 

Using various theoretical models, theorists can demonstrate quantitatively the halo den- 
sity distributions. Furthermore, the giant halo nuclei with more halo nucleons have been 
predicted in exotic Zr nuclei with A > 122 using RCHB method [11]. The prediction of giant 
halo arouses a great interest of many nuclear experimenters. Unfortunately, the exotic Zr 
nuclei are too heavy to produce in nowadays accelerator. Here in this article, we investigated 
the ground state properties of the whole Ca chain by RCHB theory, and predict that giant 
neutron halo phenomena would lie in exotic Ca nuclei with A > 60 as well as the exotic Zr 
nuclei. Compared to Zr isotopes, exotic Ca nuclei have lighter mass number and would be 
much easier to synthesize. The heaviest Ca isotope now known experimentally is °’Ca, so 5 


more neutrons are needed to form the fringe giant halo nucleus ?Ca, and 9 more neutrons 
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to form the typical giant halo nucleus Ca discussed above. 

In order to confirm whether there exists a possibility of giant halo nuclei in more wider 
mass region, in Fig. 19, we show the two-neutron separation energies Son for even-neutron 
Ne, Na, Mg, and Al nuclei in drip line region. Open symbols represent the values calculated 
from the RCHB theory with the NL-SH parameter set while corresponding solid symbols 
represent the data available. We note here that the calculations were performed by keeping 
the spherical shape and hence the comparison with the available experimental data are to 
be made with care. In Fig. 20, we plot Sən for exotic even-neutron Ar, K, Ca, Sc and 
Ti isotopes. From these two figures, the two-neutron separation energies Son for all these 
isotope chains are almost parallel. There are more than one line lying within 2 MeV in the 


_ drip line region. Therefore there are quite a large mass region where a giant halo may exist. 


= FIG. 19. The two-neutron separation energies S2, for even-neutron Ne, Na, Mg, and Al nuclei 
rr in the drip line region. Open symbols represent the values calculated with the RCHB theory with 
the NL-SH parameter set, while corresponding solid symbols represent the data available. The 


horizontal line at 2MeV denotes the upper limit for the possible occurrence of halos. 
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FIG. 20. Same as Fig.19, but for even Ar, K, Ca, Sc, and Ti nuclei in the drip line region. 


Many nuclear physicists are making much effort to search the drip line for heavier ele- 
l ~ ments. A new experiment [32,33] has been reported that both 37Na and **Ne are bound, 
B while 33Ne and *°Na are missed in the experiment. Looking at Fig. 19, the 37Na and *4Ne 


© nuclei lie in the drip line region and approach to the giant halo nuclei. Therefore we suggest 


5 that much more experimental effort should be devoted to extend the study in this area and 
= measure the mass to find the possible giant halo nuclei. Similar situation can also be seen 


x in Fig.20 for Ar, K, Ca, Sc and Ti isotopes near the neutron dripline. 


V. SUMMARY 


We have investigated the ground state properties of even-even proton magic O, Ca, Ni, 
Zr, Sn, and Pb isotopes with the relativistic continuum Hartree-Bogoliubov (RCHB) theory. 
We have found good agreement with available experimental data for the binding energies 
and the nuclear radii. We have shown the binding energies, Æp, two neutron separation 
energies, Sən, and root mean square radii. The predicted neutron drip nuclei are 7°O, Ca, 
IBNI, 6Zr, 16Sn, and 7°Pb, respectively. Particularly the giant halos in the neutron rich 
Ca and Zr isotopes close to the neutron drip line are predicted. The giant halo properties 
in exotic Ca isotopes with A > 60 is discussed in detail. 

Giant neutron halos in exotic Ca isotopes have been studied from the analysis of Sən, 


radii, nucleon density distribution, single particle energy levels, occupation probabilities and 
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the contributions from the continuum. The spin-orbit splitting and the potential diffuseness 
in Ca isotopes have also been investigated. Summarizing the present investigation, we 


conclude: 


1. Based on the analysis of two-neutron separation energies Son, rms radii, single-particle 
levels spectra, the orbital occupation, and the contribution of the continuum, giant halo 
phenomena are suggested to appear in Ca isotopes with A>60. Similar phenomena 


can also be seen for nuclei near the Na or Ar isotopes near the neutron dripline. 


2. The neutron drip line nucleus for Ca is Ca instead of Ca, which is caused by the 


disappearance of the N = 50 magic number due to the halo effect of the 3s1/2 orbit. 


ye 3. The giant halos developed in these nuclei are due to the pairing correlation and the 


contribution from the continuum e.g. the 3s1/2 orbit. 


4. The spin-orbit splitting in Ca isotopes decreases monotonically from the proton drip 


line to the neutron drip line for most cases. This tendency mainly comes from the 


diffuseness of nuclear potential with the neutron number. 


In this paper, the proton magic even-even nuclei from the proton drip line to the neu- 


` = tron drip line are studied in detail in RCHB theory. The important contribution from the 
"a , continuum due to the pairing correlations has been taken into account. The powerfulness of 


> the RCHB method has been demonstrated for the proton magic nuclei using the assumption 


il 


e of the spherical shape. One of the other important degrees of freedom for the exotic nuclei 
Q is the deformation. The theoretical framework for exotic nuclei with the deformation and 


contribution from continuum are in progress and will be completed soon. 
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